#!/usr/bin/env python
# -*- coding: iso-8859-1 -*-
from sage.all import *
from curves import *
from experimental import *

Pxy = PolynomialRing(GF(7), ['x', 'y'])
Px = PolynomialRing(GF(7), 'x')
(k, y) = Pxy.gens()
x = Px.gen()
fstr = 'x**5 -5*x**3 + 4 * x'
r = plotCurve(fstr)
f = eval(fstr)
h = 0 * x
C = HyperellipticCurve(f)
pontosCurva = C._points_cache_sqrt()
specials = specialPoints(pontosCurva, f, h)
polyOriginal = y ** 2 * (y + k)
polyReduzido = Pxy(f) * (y + x)


zero = 0 * x
one = 1 + zero
J = C.jacobian()
ponto = J([x, zero])
# searchSame(ponto, Px, J, f, h)
print C
r.save('myfile.png')
